Predicting economic conditions

ABSTRACT

Computer-implemented methods for identifying or assessing any type of risk and/or opportunity that may arise can include either, alone or in combination, band pass filtering, principal component analysis, random matrix theory analysis, synchronization analysis, and early-warning detection. Each technique can also be viewed as a process that takes a set of inputs and converts it to a set of outputs. These outputs can be used as inputs for a subsequent process or the outputs may be directly actionable for formulating certain economic predictions to make certain decisions.

This application is a continuation-in-part of co-pending U.S. patentapplication Ser. No. 14/078,197 entitled “Predicting EconomicConditions” and filed on Nov. 12, 2013, the entire disclosure of whichis hereby incorporated by reference.

BACKGROUND

Modern global finance depends on a high level of connectivity amongfinancial institutions. In stable market conditions such connectionsallow capital to flow freely with little regard to geography.Transactions among investors/savers and liquidity providers andliquidity users can be globally efficient. The same networked structuremay, however, become a channel of economic instability and distressamplification during global, financial stress. In the case of banks,where savers (deposits) fund investments (loans), a devaluation of loanscan create abnormal liquidity demands from deposits. In the same way, inthe case of shadow banks, where savers (money markets) fund investments(equity, debt and derivatives), a devaluation of collateral can alsocreate an abnormal liquidity demand. These abnormal liquidity demandscan escalate due to the interconnectedness of the economic system.

Studying systemic risk as compared with the emphasis on developingconventional risk management techniques in individual entities may behelpful in understanding economic conditions. The losses resulting fromsystemic risk may be taken into account. Focusing on the individual firmlevel in risk management may not be enough in managing the risk of acomplex, interconnected system of companies. A larger, systemicperspective may be desired. Therefore, micro-macro connections can bestudied, though individual components unique to a given microeconomicsystem. The topological properties of the networks in macroeconomicsystem share similarities with universal organizing principles.

BRIEF SUMMARY

The following presents a simplified summary of various aspects describedherein. This summary is not an extensive overview, and is not intendedto identify key or critical elements or to delineate the scope of theclaims. The following summary merely presents some concepts in asimplified form as an introductory prelude to the more detaileddescription provided below.

Computer-implemented methods for identifying or assessing any type ofrisk and/or opportunity that may arise can include either, alone or incombination, band pass filtering, principal component analysis, randommatrix theory analysis, synchronization analysis, and early-warningdetection. Each technique can also be viewed as a process that takes aset of inputs and converts it to a set of outputs. These outputs can beused as inputs for a subsequent process or the outputs may be directlyactionable for formulating certain economic predictions to makeinvestment decisions and the like.

BRIEF DESCRIPTION OF THE DRAWINGS

A more complete understanding of the present disclosure and theadvantages thereof may be acquired by referring to the followingdescription in consideration of the accompanying drawings, in which likereference numbers indicate like features, and wherein:

FIG. 1 illustrates one example of a network architecture and dataprocessing device that may be used to implement one or more illustrativeaspects discussed herein.

FIG. 2 illustrates a flow diagram for an exemplary process disclosedherein.

FIG. 3 illustrates a flow diagram for another exemplary processdisclosed herein.

FIG. 4 illustrates an exemplary diagram relating to one or more aspectsof the disclosure herein.

FIG. 5 illustrates a flow diagram for another exemplary processdisclosed herein.

FIGS. 6A-6C illustrate an exemplary data filtering technique disclosedherein.

FIG. 7 illustrates a graph relating to an exemplary process of FIG. 5.

FIG. 8A illustrates a flow diagram for another exemplary filteringprocess disclosed herein.

FIG. 8B illustrates an exemplary graph relating to the exemplary processof FIG. 8A.

FIG. 9 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 10 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 11 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 12 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 13 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 14 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 15 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 16 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 17 illustrates another exemplary graph relating to principlecomponent analysis.

FIG. 18 illustrates an exemplary flow chart pertaining to principlecomponent analysis.

FIG. 19 illustrates a flow chart pertaining to an exemplary datafiltering technique and principle component analysis.

FIG. 20A illustrates a flow chart pertaining to an exemplary randommatrix theory analysis.

FIG. 20B illustrates another flow chart pertaining to an exemplaryrandom matrix theory analysis.

FIGS. 21A-21H show the evolution of exemplary economic data of anexemplary system over time.

FIG. 22A shows the evolution of an exemplary order parameter over time.

FIG. 22B illustrates an exemplary phase distribution of economic data.

FIGS. 23A-23D illustrate exemplary Home Price Index data and derivativesthereof.

FIG. 24 shows an exemplary process for early warning detection ofeconomic conditions.

FIG. 25 illustrates an apparatus with a data processor and an earlywarning system in accordance with an aspect of the embodiments.

FIG. 26 illustrates a computer system that supports processing ofinformation about a real-world system and generating an early warningsignal in accordance with an aspect of the embodiments.

FIG. 27 illustrates a flow diagram for generating an early warningsignal from processed information and responding to the early warningsignal in accordance with an aspect of the embodiments.

FIG. 28 illustrates a flow diagram for generating an early warningsignal from processed information and responding to the early warningsignal in accordance with an aspect of the embodiments.

FIG. 29 illustrates a flow diagram for sending an early warning signalto a plurality of computing devices in accordance with an aspect of theembodiments.

FIG. 30 illustrates a flow diagram for displaying information obtainedfrom an early warning signal at a computing device in accordance with anaspect of the embodiments.

DETAILED DESCRIPTION

In the following description of the various embodiments, reference ismade to the accompanying drawings, which form a part hereof, and inwhich is shown by way of examples various embodiments in which thedisclosure may be practiced. It is to be understood that otherembodiments may be utilized and structural and functional modificationsmay be made without departing from the scope of the present disclosure.The disclosure is capable of other embodiments and of being practiced orbeing carried out in various ways. Also, it is to be understood that thephraseology and terminology used herein are for the purpose ofdescription and should not be regarded as limiting. Rather, the phrasesand terms used herein are to be given their broadest interpretation andmeaning. For example, the use of “including” and “comprising” andvariations thereof is meant to encompass the items listed thereafter andequivalents thereof as well as additional items and equivalents thereof,and the use of the terms “mounted,” “connected,” “coupled,”“positioned,” “engaged” and similar terms, is meant to include bothdirect and indirect mounting, connecting, coupling, positioning andengaging.

As noted above, various aspects of the disclosure relate to predictingeconomic conditions. Before discussing these aspects in greater detail,however, several examples of a network architecture and a dataprocessing device that may be used in implementing various aspects ofthe disclosure will first be discussed.

Detailed Description of Example Network Architecture and Data ProcessingDevice that May be Used to Implement Methods for Determining EconomicConditions:

FIG. 1 illustrates one example of a network architecture and dataprocessing device that may be used to implement one or more illustrativeaspects. Various network nodes 103, 105, 107, and 109A-F may beinterconnected via a wide area network (WAN) 101, such as the Internet.Other networks may also or alternatively be used, including privateintranets, corporate networks, LANs, wireless networks, personalnetworks (PAN), and the like. Network 101 is for illustration purposesand may be replaced with fewer or additional computer networks. A localarea network (LAN) may have one or more of any known LAN topology andmay use one or more of a variety of different protocols, such asEthernet. Devices 103, 105, 107, 109A-F and other devices (not shown)may be connected to one or more of the networks via twisted pair wires,coaxial cable, fiber optics, radio waves or other communication media.For example, the above connections can be made via the internet, bluetooth, WiFi, infrared, or any other known method of wirelesstransmission.

As shown in FIG. 1, devices 109A-F may include personal computers suchas desktops, laptops, notebooks, mobile telephones or smartphones withapplications and other functionality, a handheld device with Wi-Fi orother wireless connectivity (e.g., wireless enabled tablets, tabletcomputers, PDAs, and the like), displays with built-in or externalmemories and processors, or any other known computer, computing device,or handheld computer can also be connected to one or more of thenetworks described herein. It is also contemplated that other types ofdevices such as ATMs, kiosks, and other cash handling devices can beconnected to one or more of the networks described herein. These devicescan be enabled to communicate with wireless access points which in oneexample can be a series of cellular towers hosted by a service provider.Additionally, the wireless access points may be Wi-Fi (e.g., compatiblewith IEEE 802.11a/b/g/ and the like wireless communication standards)connections and the computing devices may obtain access to the Internetat these connections. Other known techniques may be used to allowdevices to connect with a network.

The term “network” as used herein and depicted in the drawings refersnot only to systems in which remote storage devices are coupled togethervia one or more communication paths, but also to stand-alone devicesthat may be coupled, from time to time, to such systems that havestorage capability. Consequently, the term “network” includes not only a“physical network” but also a “content network,” which is comprised ofthe data—attributable to a single entity—which resides across allphysical networks.

The components may include data server 103, web server 105, and clientcomputers 107, and devices 109A-F. Data server 103 provides overallaccess, control and administration of databases and control software forperforming one or more illustrative aspects as described herein. Dataserver 103 may be connected to web server 105 through which usersinteract with and obtain data as requested. Alternatively, data server103 may act as a web server itself and be directly connected to theInternet. Data server 103 may be connected to web server 105 through thenetwork 101 (e.g., the Internet), via direct or indirect connection, orvia some other network. Users may interact with the data server 103using remote computers 107, devices 109 A-F, e.g., using a web browserto connect to the data server 103 via one or more externally exposed websites hosted by web server 105. Client computers 107, 109 may be used inconcert with data server 103 to access data stored therein, or may beused for other purposes. For example, from client device 107 or devices109 A-F a user may access web server 105 using an Internet browser, asis known in the art, or by executing a software application or app thatcommunicates with web server 105 and/or data server 103 over a computernetwork (such as the Internet).

Servers and applications may be combined on the same physical machines,and retain separate virtual or logical addresses, or may reside onseparate physical machines. FIG. 1 illustrates just one example of anetwork architecture that may be used, and those of skill in the artwill appreciate that the specific network architecture and dataprocessing devices used may vary, and are secondary to the functionalitythat they provide, as further described herein. For example, servicesprovided by web server 105 and data server 103 may be combined on asingle server.

Each component 103, 105, 107, 109 may be any type of known computer,server, or data processing device as discussed herein. Data server 103,e.g., may include a processor 111 controlling overall operation of thedata server 103. Data server 103 may further include RAM 113, ROM 115,network interface 117, input/output interfaces 119 (e.g., keyboard,mouse, display, printer, or the like.), and memory 121. I/O 119 mayinclude a variety of interface units and drives for reading, writing,displaying, and/or printing data or files. Memory 121 may further storeoperating system software 123 for controlling overall operation of thedata processing device 103, control logic 125 for instructing dataserver 103 to perform aspects as described herein, and other applicationsoftware 127 providing secondary, support, and/or other functionalitywhich may or may not be used in conjunction with one or more aspectsdescribed herein. The control logic may also be referred to herein asthe data server software 125. Functionality of the data server softwaremay refer to operations or decisions made automatically based on rulescoded into the control logic, made manually by a user providing inputinto the system, and/or a combination of automatic processing based onuser input (e.g., queries, data updates, or the like).

Memory 121 may also store data used in performance of one or moreaspects, including a first database 129 and a second database 131. Insome embodiments, the first database may include the second database(e.g., as a separate table, report, or the like). That is, theinformation can be stored in a single database, or separated intodifferent logical, virtual, or physical databases, depending on systemdesign. Devices 105, 107, 109 may have similar or different architectureas described with respect to device 103. Those of skill in the art willappreciate that the functionality of data processing device 103 (ordevice 105, 107, 109 A-F) as described herein may be spread acrossmultiple data processing devices, for example, to distribute processingload across multiple computers, to segregate transactions based ongeographic location, user access level, quality of service (QoS), or thelike.

One or more aspects may be embodied in computer-usable or readable dataand/or computer-executable instructions, such as in one or more programmodules, executed by one or more computers or other devices as describedherein. Generally, program modules include routines, programs, objects,components, data structures, or the like that perform particular tasksor implement particular abstract data types when executed by a processorin a computer or other device. The modules may be written in a sourcecode programming language that is subsequently compiled for execution,or may be written in a scripting language such as (but not limited to)HTML or XML. The computer executable instructions may be stored on acomputer readable medium such as a hard disk, optical disk, removablestorage media, solid state memory, RAM, or the like. As will beappreciated by one of skill in the art, the functionality of the programmodules may be combined or distributed as desired in variousembodiments. In addition, the functionality may be embodied in whole orin part in firmware or hardware equivalents such as integrated circuits,field programmable gate arrays (FPGA), and the like. Particular datastructures may be used to more effectively implement one or moreaspects, and such data structures are contemplated within the scope ofcomputer executable instructions and computer-usable data describedherein.

Detailed Description of Example Methods and Systems for DeterminingEconomic Conditions:

Identifying structural attributes of the interconnected world economicnetwork can provide clues about which characteristics of the systemcorrelate with a certain degree under systemic distress leading toeconomic instability or economic stability. Based on the understandingderived from these studies, an early-warning system ofeconomic-instability detection can be developed. This early-warningsystem can also be used in identifying and assessing any type of riskand opportunity that may arise. It can be used to mitigate risk ormaximize gain.

Additionally, aggregate macroeconomic events at the largest scale arebuilt up from individual microeconomic agents or activities at the mostgranular level. This may be reflected in utilizing network analysis toidentify economic instability in financial systems.

The interlinked ebb-and-flows of financial markets, which show thedynamic structure of financial systems, affect systemic risk inherent init. In an ideal situation where relevant data or their possible proxiesare available, banking systems that are composed of a number of severalconnected banks can be constructed. Key parameters that define thestructure of the system can be varied accordingly. These may include thelevel of capitalization, the degree to which banks are connected, thesize of interbank exposures and the degree of concentration of thesystem. The result of the analysis of the influence of these parameterscan be used to potentially modify the financial system and mitigaterisks as needed.

Network analysis can be applied to analyze the increasingly complex andglobally interlinked financial markets. A network is simply a collectionof points (or nodes) joined by lines (or edges). Networks provide asimple but useful representation of real-world systems of interactingcomponents. The internet, for instance, can be represented as a networkof computers linked by data connections. Other examples include socialnetworks of friendships between individuals and networks of businesscontacts. Networks are increasingly common in the study of biology,epidemiology, genetics, transportation, computer software, and so on.Networks can also be applied in exploring some possible applications ofnetwork analysis in financial systems.

It is possible to model the propagation of failures in a financialsystem as an epidemic spreading process in a network of interlinkedstocks and flows of money. Starting, for example, from a small number offailed banks, the aim is to characterize the probability that failurespropagate at the systemic level as a function of some relevantparameters, like the connectivity of the network and concentrations,correlations, sensitivities, leverage and liquidity of asset andliability classes at banks. One of the key problems in this context isto understand the role of the network structure in relation to economicinstability.

System variables that show purely random fluctuations or perfectlyperiodic rhythms define idealized extremes. In fact, some parameterswhose fluctuations may seem random when viewed in isolation can behighly predictable in the temporal context of variation in otherparameters; hence, variations in one system variable can conveysubstantial information about variation in another. An organism (or acompany) that is capable of learning these correlations can exploit themin order to anticipate vital changes in the system. Stressful stimulimay be important not because of their immediate and direct consequences,but in the information they convey about the overall state of the systemand its likely trajectory. Therefore, survival may depend much on theart of correctly identifying the system's trajectory coupled with theability to adapt to changes in the system.

Robustness of the economic system and the ability to withstand economicshock can be promoted by modularity—the degree to which nodes of asystem can be decoupled into relatively discrete components. Forexample, a basic principle in management of forest fires and epidemicsis that if there is strong interconnection among the system elements, aperturbation will encounter nothing to stop it from spreading. But ifthe system is appropriately compartmentalized through, say, introductionof fire breaks and quarantining, disturbance or risk is more easilycountered.

The introduction of modularity, however, will often involve a trade-offbetween local and systemic risk. Moreover, the wrongcompartmentalization in financial markets could preclude stabilizingfeedbacks, such as mechanisms for maintaining liquidity of cash flowsthrough the financial system, where fragmentation leading to illiquiditycould actually increase systemic risk. Redundancy of components andpathways, in which one can substitute for another, is also a key elementin the robustness of complex systems, and effective redundancy is notindependent of modularity.

There are only coarse or indirect options for control of the financialsystem. The tools available to policymakers are designed to modifyindividual incentives and individual behaviors in ways that will supportthe collective good. Such top-down efforts to influence individualbehaviors can often be effective, but in certain instances it may be achallenge to control the spread of synchronized behaviors or to managefinancial crises in an optimal way.

Although the study of payment flows is of interest to central bankers,in certain instances it may miss aspects of systemic risk, namely publicperception and asset valuation associated with the interaction ofcounterparties (the mutual financial obligations and exposures that linkcompanies). Such company networks are helpful in studying the effects ofinflated asset-pricing, credit crises and the poorly understood butpotentially worrying effects of the current widespread use ofderivatives by investment banks to manage risk in real-time. Whateverthe case, it seems that networks that define financial reality andglobal markets may be of use to understanding the market robustness andability to thrive as well as its potential vulnerability to collapse.

Complex dynamical systems, financial system included, may havetransitional points where a sudden shift to a totally differentdynamical regime may take place. Though predicting such tipping pointsis very difficult, the generic early-warning signs that may indicate ifa critical threshold is approaching may exist in reality.

Many complex dynamical systems have critical thresholds called tippingpoints where the system abruptly shifts from one state to another. Inglobal finance, there is big concern about a sudden systemic decline inmarket prices that can threaten the global financial system which, inturn, may lead to global economic downturns. It is traditionallydifficult to predict such critical transitions since the state of thesystem may show little change before the tipping point is reached.However, it seems that certain generic signs may take place in a wideclass of systems as they approach a critical point. The dynamics ofsystems near a critical point have generic properties, regardless ofdifferences in the details of each system. Critical thresholds for suchtransitions correspond to bifurcations. In the so-called catastrophicbifurcation, once a threshold is exceeded, a positive feedback mechanismpushes the system through a phase of directional change towards acontrasting state. The transitions from a stable equilibrium to a cyclicor chaotic state can also happen in some other types of bifurcations.

Three possible early-warnings in the dynamics of a system approaching atipping point may be used in emerging risk detection: slower recoveryfrom perturbations, increased autocorrelation and increased variance.

Because slowing down causes the intrinsic rates of change in the systemto decrease, the state of the system at any given moment becomes moreand more like its past state. The resulting increase in “memory” of thesystem can be measured in various ways from the frequency spectrum ofthe system. One approach is to look at the lagging autocorrelation,which can directly be interpreted as slowness of recovery in suchnatural perturbation regimes. Analysis in the models exposed tostochastic forcing confirms that if the system is driven graduallycloser to a catastrophic state, there is a marked increase inautocorrelation that builds up long before the critical transition takesplace. This is true even for realistic models. Increased variance in thepattern of fluctuations can be seen as another possible consequence ofcritical slowing as a critical transition is approached. As theEigenvalue approaches zero, the impacts of shocks do not decay, andtheir accumulating effect increases the variance of the state variable.In principle, critical slowing could reduce the ability of the system totrack the fluctuations, and thereby produce an opposite effect onvariance. However, analysis shows that an increase in the varianceusually arises and may be detected before a critical transition takesplace. This detection technique can be implemented both temporally andspatially.

Methods for identifying or assessing any type of risk and/or opportunitythat may arise can include either, alone or in combination, band passfiltering, principal component analysis, random matrix theory analysis,synchronization analysis, and early-warning detection. Each techniquecan also be viewed as a process that takes a set of inputs and convertsit to a set of outputs. These outputs can be used as inputs for asubsequent process or the outputs may be directly actionable forformulating certain economic predictions to make investment decisionsand the like.

A general diagram of this approach is illustrated in FIG. 2. FIG. 2shows a diagram of an analysis where the system analyzes a series ofinputs, which can be any financial data, for example, consumerleverages, commercial leverages, and unemployment rates and the like.These inputs are processed using any of the techniques described hereinto produce an output or economic indicator. These economic indicatorscan be used to develop strategies for making investments, mitigatingrisks, creating new products, and the like. The various outputs from thetechniques described herein can be used to arrive at a decision, e.g. tosell, invest, mitigate risk, and the like.

As schematically illustrated in FIG. 3, either approach works byconverting various sets of inputs into outputs that may be used to guidefinancial decisions. Based on the understanding derived from thesestudies, an early-warning system of assessing and dealing with variousrisks (including, but not limited to, economic instability) is alsodeveloped. In one example, as shown in FIG. 3, various outputs can bereviewed to arrive at the financial decision, for example, outputs froma top-down review of financial data, bottom up review of financial data,and other outputs.

Certain techniques can be used in a top-down approach that reviewsinputs from the macro level down to the micro level or a bottom upapproach that works in an opposite fashion. Considerable noises makemodeling time series data in its most granular level (such as modelingaccount level in portfolio analysis or relationship networks in finance)more difficult. To deal with these difficulties, a top-down approach canbe conducted. This approach is intended to firstly understand thebehavior at the most aggregate level (top), where the relevant nationallevel macroeconomic time series are used. Then the analysis is graduallyperformed towards the more granular level (down), where time series databeing used includes, but not limited to, relevant geographic or statelevel, and segmental level. As information is viewed down from the mostaggregate level to the more granular level, any useful derivableinformation is taken into consideration, such as how shocks, policies,and strategies may contribute to the behaviors being observed.

A complexity pyramid composed of various levels of granularity is shownin FIG. 4. The bottom of the pyramid shows the representation ofspecific individual to small-sized economic agent network or circuit(level 1). Insights into the logic of most granular organization can beachieved when the economic system is viewed as a complex network inwhich the components are connected by some kind of relationship links.At the lowest level, these components form medium-sized networks (level2), which in turn are the building blocks of functional modules (level3). These modules are nested, generating scale-free or other kinds ofnetwork architectures (level 4). Although individual components areunique to a given microeconomic system, the topological properties ofthe networks share similarities with universal organizing principles.

Band Pass Filtering Technique:

A band pass filtering technique can be applied to economic inputs inorder to extract a signal from noise by removing certain data orfrequencies that are too rapid or too slow to be “true,” similar toelectronic signal processing. Certain filters can be applied to economictime series data to identify the particular phase in an economic cycle.The analysis may help guide certain decisions such as whether to loosenor tighten credit availability.

FIG. 5 illustrates a process flow of an application of a band passfilter to a series of inputs 502. These inputs can be consumer andcommercial leverages and unemployment rates, and can be eithermacroeconomic or microeconomic data. First the series of inputs iscollected and transformed to a proper frequency domain at 504. The bandpass filter can then be applied at step 506 to generate a meaningfulcycle signal in step 508. FIGS. 6A-6C show a raw time series transformedinto a frequency domain with the signal extracted. In particular, a bandpass filter is applied to the data depicted in FIG. 6A and FIG. 6B toproduce the filtered data shown in FIG. 6C. The data shown in FIG. 6 canthen be evaluated to determine a synchronization point or region of thedata that can be used to make economic predictions based on thesynchronization events.

As shown in FIG. 7, the band pass filtering technique can be applied tomultiple economic indicators simultaneously. If the correct band passfilter is implemented, when each of the signals synchronize, it canindicate a transition from an equilibrium state to either a positivestate (economic boom) or a negative state (recession). FIG. 7illustrates a synchronization event in the period around 2010, whichshows a positive state of the economy.

During a synchronization event, the risk or reward can increase and thebenefits of diversification are diminished. The synchronization levelcan be measured based on the collective correlation of a group of timeseries. This correlation could also be applied to additional measuresother than the simple filtered series.

For example, the absolute value of the volatility of the changes can bedetermined. This process is depicted generally in FIG. 8A. As shown inFIG. 8A, time series data 802 can be processed by a band pass filter, atstep 804, can be applied to the economic data and the absolute value ofthe data can be taken at step 806 to locate correlations across all ofthe series of data at step 808. The correlation and concentrationmeasures can then be reviewed at step 810 to predict certain economictrends, for example, the level of leverage of consumer and commercialsectors. An example output of the absolute values is illustrated in FIG.8B. As shown in 8B, in this example, a higher correlation value may beassociated with a weaker state in the economy, whereas a lowercorrelation value may indicate a stronger performing economy. However,higher correlation value may reflect a better state of the economy, ifthe level of saving is considered. Therefore, higher correlation may beassociated with stronger or weaker performing economy based on the databeing observed.

Principal Component Analysis:

Once it is determined whether or not the economy is moving in a positiveor negative direction, the next step is to determine if there are anyparticular elements driving the change. On a macro level this could becertain industries or asset classes. On a micro level it could be acertain company's economic performance or group of companies' economicperformances. Principal component analysis or PCA can be used toidentify the particular “drivers” of the trends. This can indicate whichsubcomponents are accounting for the change in state.

As shown in FIG. 9, largest ratios of two predetermined components mayillustrate a particular economic driver. As shown in FIG. 9, twosubcomponents can be compared by calculating the ratio between theirfirst and second principal components. In one example, the firstprincipal component can represent systemic risk, for example, arecession or other event external to the data being studied. The secondprincipal component can be non-systemic (idiosyncratic) risk, which canbe specific to a particular subcomponent. In FIG. 9, Subcomponent 5 maybe the market driver because the ratio of the first principal componentto the second principal component is the largest.

In another example, certain economic instability markers can beidentified using spectral analysis and principal component analysis. Inthis case, because the typical size of a U.S. business cycle isapproximately 70 months, this data can be analyzed as part of a study.Spectral analysis of correlation matrices can be conducted to determineeconomic instability markers.

In certain examples, the analysis can be applied to uncollectables,where it has been declared that certain amounts of money are unlikely tobe collected, by sectors as shown in FIGS. 10-13. This analysis can alsobe applied to the Home Price Index (“HPI”) by states as shown in FIGS.14-17. By first scaling the time series, a correlation matrix isgenerated within a moving time window with the length of roughly onebusiness cycle (˜70 months). In each of the moving windows, principalcomponent analysis (PCA) can be conducted. The first principal component(PC1) is usually called “systemic risk,” which basically represents thecollective response of all component time series under consideration toexogenous stimuli. The behavior of the PC1 eigenvector contents (orfactor loadings) over time can be investigated. A quality calleddispersion level is provided to give a measure of disparity or diversityof the time series being analyzed.

For each analysis, a set of figures are shown: the time series PCA(FIGS. 10 and 14), PC1 Eigenvector contents (FIGS. 12 and 16), and PC1dispersions (FIGS. 13 and 17). The vertical black ribbons showrecessions, and the vertical green ribbons show the 1995/1996 periodwhere the bad consumer cycle took place without being accompanied by arecession. As can be seen in FIGS. 10 and 14, the time series data canbe scaled. The eigenvalues obtained using PCA are shown in FIGS. 11 and15 and can be normalized such that their total summation is equal to100%. FIGS. 12 and 16 show the contents of the first principal component(or, principal component number 1).

Exemplary charts of the dispersion level over time are shown in FIGS. 13and 17. The dispersion level may be taken as representing the level ofeconomic instability risk of the system under consideration. From thedynamical behavior of the dispersion level, an early-warning tool ofeconomic instability risk may be developed. In particular, the higherthe value of the dispersion level the less contagious the system is. Asshown in FIGS. 13 and 17 recessions promoted the decrease of thedispersion level. This can be interpreted as the phenomenon of theeconomic instability of systemic risk.

FIG. 18 shows an exemplary flow chart for implementing an algorithm thatmay detect economic instability risk based on eigenvector dispersion“level.” In step 1800, principal component analysis is applied toseveral time series under consideration by using an appropriate timewindow. In step 1802, the eigenvector content of principal componentnumber 1 is normalized for each time window. The eigenvectors are thensorted, e.g., from smallest to largest, say from 1 to 8 for each timewindow in step 1804. In step 1806, the distances between theeigenvectors are then calculated e.g. 1 to 2, 2 to 3, 3 to 4, and so on,specific for each time window (note that the order of eigenvectors maychange for different time window). In step 1808, for each time window,the mean (or, the variance) of the remaining distances is calculated,where this entity represents the dispersion “level” at specific timewindow (the higher the value, the higher the dispersion level). In step1810, economic instability is then predicted based on the dispersionlevel.

Band Pass Filtering and Principal Component Analysis:

Band pass filtering and principal component analysis can be used inconjunction with each other to make economic predictions. FIG. 19 showsan exemplary flow chart of an example implementation of band passfiltering and principal component analysis being used together todetermine economic instability risk. In this example, publicly availablemacroeconomic time series data 1902 can be applied to a band pass filter1904 to obtain a filtered data series 1906. The filtered data series canthen be analyzed to determine a particular economic cycle state, e.g.positive or negative trend at step 1908. Once a state of the market isdetected (steps 1910-1914), principle component analysis can be appliedto determine the particular driver of the state of the market at step1916. Using the drivers obtained in the analysis (as well as thecorrelation/concentration) at step 1918, a portfolio optimizer 1922 can,for example, be implemented for making investments, mitigating risks,creating new products, and the like at step 1924. Additionally,proprietary micro level data can be analyzed at step 1920.

Random Matrix Theory:

Another approach to spectral analysis of economic data is using RandomMatrix Theory (RMT) to analyze economic data and predict trends in themarket. Random Matrix Theory (RMT) was initially proposed by Wigner andDyson in 1960s for studying the spectrum of complex nuclei. RMT can beused to identify and model phase transitions and dynamics in physicalsystems and can also be used to create financial and economical models.

For example, RMT can be used to estimate the number of dimensions(components) of a data correlation matrix by comparing the statistics ofthe observed eigenvalues, i.e. the eigenvalues of the data correlationmatrix to those of a random matrix counterpart. The density distributionof eigenvalues for such random matrices is known, so that the comparisonbetween the observed eigenvalues and the analytical “null” distributionof RMT can be used to obtain an estimate of the number of components.Specifically, the number of observed eigenvalues larger than theanalytical maximum provides an estimate of the number of significanteigenvalue components. Hence, the observed eigenvalues of the datacorrelation matrix larger than the theoretical maximum provides areasonable approximation to the number of principal components toretain, i.e., the number of eigenvalue components obtained using PCAthat really matter in the calculation. RMT is analytical in nature, andthe computational cost of RMT may be small in certain instances.

Taking a matrix whose elements are the correlation coefficient values ofmultiple time series of interest, the standard Pearson correlationcoefficients can be defined as:

${c( {g_{i},g_{j}} )} = {\frac{1}{N}{\sum\limits_{{k = 1},N}^{\;}\;{( \frac{g_{ik} - M_{gi}}{\sigma_{gi}} )( \frac{g_{jk} - M_{gj}}{\sigma_{gj}} )}}}$where M_(gi) and M_(gj) are the average of g_(i) and g_(j) respectively,and σ_(gi) and σ_(gj) are their corresponding standard deviations, and Nis the total number of observations.

The statistical properties of the eigenvalues of random matrices areknown in the limit of very large dimensions. Particularly, in the limitN→∞, L→∞, such that Q≡L/N is fixed, the distribution P_(RM) (λ) ofeigenvalues λ of the random correlation matrix is given by

${P_{RM}(\lambda)} = {\frac{Q}{2\;\pi} + \frac{\sqrt{( {\lambda_{+} - \lambda} )( {\lambda - \lambda_{-}} )}}{\lambda}}$for λ within the bounds λ≦λ_(i)≦λ₊, where λ⁻ and λ₊, are the minimum andmaximum eigenvalues of the random correlation matrix, respectively,given by

$\lambda_{\pm} = {1 + {\frac{1}{Q} \pm {2{\sqrt{\frac{1}{Q}}.}}}}$

RMT focuses on the study of statistical properties of eigenvalue spacingbetween consecutive eigenvalues. From RMT, distribution of eigenvaluespacing of real and symmetrical random matrices follows two universallaws depending on the correlativity of eigenvalues. Strong correlationof eigenvalues leads to statistics described by the Gaussian OrthogonalEnsemble (GOE). On the other hand, eigenvalue spacing distributionfollows Poisson statistics if there is no correlation betweeneigenvalues. To express it differently, eigenvalue spacing distributionof a random matrix with non-zero values only for its diagonal (orblock-diagonal parts) follow Poisson statistics, because eigenvalues ofthis system are not correlated due to the absence of interaction betweendiagonal (or block-diagonal) parts.

To validate RMT estimate and to ensure that the theoretical nulldistribution does not deviate significantly from that of empirical null(deviations might be expected because of the finite size of the matrixand because the data may not be gaussian), the data matrix may bescrambled (for each row, a distinct permutation of the columns isperformed) and verified that RMT predicted zero significant components.

The nearest neighbor spacing distribution (NNSD) of eigenvalues, P(s),of RMT can also be used in some calculations, as discussed below. Thisis defined as the probability density of the so-called unfoldedeigenvalue spacing s=e_(i+1)−e_(i) where e_(i)=N_(av)(E_(i)), andE_(i)(i=1, . . . , N) is the eigenvalues of the matrix (N being theorder of the matrix), and N_(av) is the smoothed integrated density ofeigenvalues obtained by fitting the original integrated density to acubic spline or by local density average. From RMT, P(s) of the GOEstatistics closely follows Wigner-Dyson distribution P_(GOE)(s)≈½ πsexp(−πs²/4). In the case of Poisson statistics, P(s) is given by Poissondistribution P_(Poisson)(s)=exp(−s). The difference between Wigner-Dysonand Poisson distributions manifests in the regime of small s, whereP_(GOE) (s→0)=0 and P_(Poisson) (s→0)=1.

Note that RMT technique can be applied to both higher and lower levelviews. In higher level, it may be used to extract the importantcompanies/groups which are most responsible for the economic dynamicsbeing observed. At a lower level, it may be used to extract theunderlying skeleton of some complex networks.

To test the “modularity” of the clustering, lower values of correlationcoefficients can be removed as given by the equation

${c( {g_{i},g_{j}} )} = {\frac{1}{N}{\sum\limits_{{k = 1},N}^{\;}{( \frac{g_{ik} - M_{gi}}{\sigma_{gi}} )( \frac{g_{jk} - M_{gj}}{\sigma_{gj}} )}}}$from the data of interest, starting from the lowest. Using chi-squaretest, a sharp transition from a Wigner-Dyson distribution to a Poissondistribution would be observed at a certain defined “cutoff” level q.Once this takes place, the desired clustering can be obtained. It shouldbe noted that this approach is different from existing clusteringmethods, where here cutoffs or thresholds used for clustering aredetermined self-consistently by the transition given by RMT. Forsegmentation, the time series that are used could be FICO scores,geographical performance, loss, revenue, and so on.

FIG. 20A shows an example flow chart using the RMT method to determineeconomic instability in financial systems in segmentation or clustering,for example in risk segmentation or population clustering. First data ofmultiple time series of interest is obtained in step 2000. Next,standard Pearson correlation coefficients using the data are built atstep 2002. The matrix whose elements are the correlation coefficientvalues (call the matrix CCM or Correlation Coefficient Matrix) is thenformulated at step 2004. The nearest neighbor spacing distribution bygenerating it from the differences between nearest neighbor eigen-valuesderived from CCM is then obtained at step 2006. The lowest value ofcorrelation coefficients is iteratively set to zero, and the calculationis repeated at step 2008. Using a chi-square test, a sharp phasetransition from Gaussian Orthogonal Ensemble (RMT) distribution toPoisson distribution would be observed at a certain “threshold” level atstep 2010. Once this takes place, optimal segmentation or clustering isachieved. The mathematical philosophy is basically trying to make CCM asdiagonal as possible, but not necessarily completely diagonal, and theprocess should be stopped once a phase transition takes place at step2012. The output could be in form of optimal segmentation or clustering,and software could be built to detect the threshold where the transitiontakes place. This could be associated with economic instability risks orsimply optimal risk segmentation or population clustering only.

An example flow chart using RMT to estimating dimensionality of the dataused for PCA is depicted in relation to FIG. 20B. First, the data ofmultiple time series of interest is determined at step 2020. Next, thestandard Pearson correlation coefficients are formulated using the dataat step 2022. The matrix whose elements are the correlation coefficientvalues (call the matrix CCM or Correlation Coefficient Matrix) is thendetermined at step 2024. The relevant RMT is formed where the densitydistribution of eigen-values for such random matrices is known at step2026. The statistics of the observed eigen-values, i.e. the eigen-valuesof the data correlation matrix is then compared to those of the RMTcounterpart at step 2028. At step 2030, the number of observedeigen-values larger than the analytical maximum from RMT provides anestimate of the number of significant eigen-value components of thedata. The observed eigen-values of the data correlation matrix largerthan the theoretical maximum from RMT provides a reasonableapproximation to the number of principal components to retain, i.e. thenumber of eigen-value components obtained using Principal ComponentAnalysis (PCA) that are relevant in the calculation. The retainedprincipal components (whose numbers are usually much lower than theoriginal dimensionality of the data) can be used in the analysis.

Early Warning Detection of Economic Conditions:

Complex dynamical systems, financial system included, can havetransitional points where a sudden shift to a totally differentdynamical regime may take place. Though predicting such tipping pointscan be difficult, generic early-warning signs may indicate if a certainthreshold is approaching.

Many complex dynamical systems have critical thresholds called tippingpoints where the system abruptly shifts from one state to another. Forexample, in global finance, there is big concern about a sudden systemicdecline in market prices that can threaten the global financial systemwhich, in turn, may lead to a global financial crisis. It may bedifficult to predict such critical transitions since the state of thesystem may show little changes before the tipping point is reached.Additionally, models of complex systems may not be accurate forpredicting where critical thresholds are located. However, it seems thatcertain generic signs may take place in a wide class of systems as theyapproach a critical point. The dynamics of systems near a critical pointhave generic properties, regardless of differences in the details ofeach system. Critical thresholds for such transitions correspond tobifurcations. In the so-called catastrophic bifurcation, once athreshold is exceeded, a positive feedback mechanism pushes the systemthrough a phase of directional change towards a contrasting state. Thetransitions from a stable equilibrium to a cyclic or chaotic state canalso happen in some other types of bifurcations.

FIGS. 21A-21H illustrate eight panels of the same exemplary economicdata, which show the evolution of an exemplary system over time. Eachpanel in FIGS. 21A-21H show three graphs: the left graph shows afrequency distribution, the middle graph shows a phase distribution, andthe right graph shows the phases in polar form on a unit circle. Asshown in FIGS. 21A-21C, from time=0 to 20 months, the exemplary economicdata is out of phase with each other. However, starting at time=30months, the exemplary economic data starts to move into phase with eachother.

FIG. 22A shows the evolution of the order parameter over time, and FIG.22B illustrates how the phase distribution is locked to have the valueswhich are confined within a certain (moving) range only. The orderparameter, in this case, shows the level of synchronicity of the systemas a whole.

In one example, indicators of whether a system is getting close to acritical threshold may be related to a phenomenon called criticalslowing down, for example, fold catastrophe. At a fold bifurcationpoint, the dominant eigenvalue characterizing the rates of change aroundthe equilibrium is zero. Therefore, as the system approaches suchcritical points, it becomes increasingly slow to recover from smallperturbations. Such slowing down typically starts far from thebifurcation point, and that recovery rates decrease smoothly to zero asthe critical points is approached.

The recovery rate after a small perturbation will be reduced, and willapproach zero when a system moves towards a catastrophic bifurcationpoint, which can be further explained by the following simple dynamicalsystem, where y is a positive scaling factor and a and b are parameters:

$\frac{\mathbb{d}x}{\mathbb{d}t} = {{\gamma( {x - a} )}{( {x - b} ).}}$

This model has two equilibria, x ₁=a and x ₂=b, of which one is stableand the other is unstable. If the value of a equals to that of b, theequilibria collide and exchange stability (in a transcriticalbifurcation). Assuming that x ₁ is the stable equilibrium, it can bedetermined what happens if the state of the equilibrium is slightlyperturbed

${( {x = {{\overset{\_}{x}}_{1} + ɛ}} ):\frac{\mathbb{d}( {{\overset{\_}{x}}_{1} + ɛ} )}{\mathbb{d}t}} = {{f( {{\overset{\_}{x}}_{1} + ɛ} )}.}$Here f(x) is the right hand side of the above equation

$\frac{\mathbb{d}x}{\mathbb{d}t} = {{\gamma( {x - a} )}{( {x - b} ).}}$Linearizing the equation using a first-order Taylor expansion yieldswhich simplifies to

$\begin{matrix}{\frac{\mathbb{d}( {{\overset{\_}{x}}_{1} + \varepsilon} )}{\mathbb{d}t} = {f( {{\overset{\_}{x}}_{1} + ɛ} )}} \\{{{\approx {{f( {\overset{\_}{x}}_{1} )} + \frac{\partial f}{\partial x}}}❘{{\overset{\_}{x}}_{1}ɛ}},{{f( {\overset{\_}{x}}_{1} )} + \frac{\mathbb{d}ɛ}{\mathbb{d}t}}} \\{= {{{f( {\overset{\_}{x}}_{1} )} + \frac{\partial f}{\partial x}}❘ {{\overset{\_}{x}}_{1}ɛ}\Rightarrow\frac{\mathbb{d}e}{\mathbb{d}t} }} \\{= {\lambda_{1}{ɛ.}}}\end{matrix}$

With eigenvalues λ₁ and λ₂ in this case,

${\lambda_{1} = {{\frac{\partial f}{\partial x}❘a} = {- {\gamma( {b - a} )}}}},$and, for the other equilibrium

$\lambda_{2} = {{\frac{\partial f}{\partial x}❘b} = {- {\gamma( {b - a} )}}}$are obtained.

If b>a then the first equilibrium has a negative eigenvalue, λ₁, and,therefore, it is stable (as the perturbation goes exponentially to zero;

$ {{{f( {\overset{\_}{x}}_{1} )} + \frac{\mathbb{d}ɛ}{\mathbb{d}t}} = {{{{f( {\overset{\_}{x}}_{1} )} + \frac{\partial f}{\partial x}}❘ {{\overset{\_}{x}}_{1}ɛ}\Rightarrow\frac{\mathbb{d}e}{\mathbb{d}t} } = {\lambda_{1}ɛ}}} ).$It is easy to see from the above equations

${\lambda_{1} = {{\frac{\partial f}{\partial x}❘a} = {- {\gamma( {b - a} )}}}},{and}$$\lambda_{2} = {{\frac{\partial f}{\partial x}❘b} = {- {{\gamma( {b - a} )}.}}}$that at the bifurcation (b=a) the recovery rates λ₁ and λ₂ are both zeroand perturbations will not recover. Farther away from the bifurcation,the recovery rate in this model is linearly dependent on the size of thebasin of attraction (b−a). For more realistic models, this is notnecessarily true but the relation is still monotonic and is often nearlylinear.

The most direct implication of critical slowing down is that therecovery rate after small perturbation can be used as an indicator ofhow close a system is to a bifurcation point. For most natural systems,it would be impossible to monitor them by systematically observingrecovery rates. However, it can be shown that as a bifurcation isapproached in such a system, certain characteristic changes in thepattern of fluctuations are expected to take place. An importantprediction is that the slowing down should lead to an increase inautocorrelation in the resulting pattern of fluctuations.

Critical slowing down will tend to lead to an increase in theautocorrelation and variance of the fluctuations in a stochasticallyforced system approaching a bifurcation at a threshold value of acontrol parameter. The example given here illustrates why this is thecase. In certain instances it may be assumed, there is a repeateddisturbance of the state variable after each period Δt (that is,additive noise). Between disturbances, the return to equilibrium isapproximately exponential with a certain recovery speed, A. In a simpleautoregressive model this can be described with the following equations:x_(n+1)=exp(λΔt)(x_(n)−x)+σε_(n), and y_(n+1)=exp(λΔt)y_(n)++σε_(n).Here y_(n) is the deviation of the state variable x from theequilibrium, ε_(n) is a random number from a standard normaldistribution and σ is the standard deviation. If λ and Δt areindependent of y_(n), this model can also be written as a first-orderautoregressive (AR(1)) process: y_(n+1)=αy_(n)+σε_(n). Theautocorrelation α≡(λΔt) is zero for white noise and close to one for red(autocorrelated) noise. The expectation of an AR(1) process:y_(n+1)=c+αy_(n)+σε_(n) is

$\begin{matrix}{{E( y_{n + 1} )} =  {{E(c)} + {\alpha\;{E( y_{n} )}} + {E( {\sigma\; ɛ_{n}} )}}\Rightarrow\mu } \\{=  {c + {\alpha\;\mu} + 0}\Rightarrow\mu } \\{= {\frac{c}{1 - \alpha}.}}\end{matrix}$

For c=0, the mean equals zero and the variance is found to be

${{Var}( y_{n + 1} )} = {{{E( y_{n}^{2} )} - \mu^{2}} = {\frac{\sigma^{2}}{1 - \alpha^{2}}.}}$

Close to the critical point, the return speed to equilibrium decreases,implying that λ approaches zero and the autocorrelation a tends to one.Therefore, the variance tends to infinity. These early-warning signalsare the result of critical slowing down near the threshold value of thecontrol parameter.

Slowing down may cause the intrinsic rates of change in the system todecrease, and the state of the system at any given moment may becomemore and more like its past state. The resulting increase in “memory” ofthe system can be measured in various ways from the frequency spectrumof the system. One example approach is to calculate a lag-1autocorrelation, which can be directly interpreted as slowness ofrecovery in such natural perturbation regimes.

Additionally, analysis in the models exposed to stochastic forcingconfirms that if the system is driven gradually closer to a catastrophicbifurcation, there is a marked increase in autocorrelation that buildsup long before the critical transition takes place. This is true evenfor realistic models. Whereas increased variance in the pattern offluctuations can be seen as another possible consequence of criticalslowing down as a critical transition is approached. As the eigenvalueapproaches zero, the impacts of shocks do not decay, and theiraccumulating effect increases the variance of the state variable. Inprinciple, critical slowing down could reduce the ability of the systemto track the fluctuations, and thereby produce an opposite effect onvariance. However, analysis shows that an increase in the varianceusually arises and may be detected before a critical transition takesplace.

Phenomenon of critical slowing down may lead to three possibleearly-warning signals in the dynamics of a system approaching abifurcation that may be used in emerging risk detection: slower recoveryfrom perturbations, increased autocorrelation, and increased variance.

An exemplary process is depicted in FIG. 24, which shows an exemplaryprocess for determining an early warning signal for economicinstability. In the first step, a data series of interest is obtained instep 2400 to determine the particular early warning signal. Next, theAuto Regressive (“AR(1)”) Model coefficients of a time series can becalculated in step 2402. The scaled lag-1 auto-correlation of thecoefficients can be determined in step 2404. Finally, the scaledsmoothed derivative of the previous computation can be determined atstep 2406. The resulting derivative can be used as an early warningdetector. For example, an indicator as a positive or negative warningcould be provided to the user. In particular, the system could beconfigured to output a positive indicator when the scaled smoothedderivative is above a predetermined threshold value and can beconfigured to output a negative indicator when the scaled smoothedderivative is below a predetermined threshold value.

In one example, the above-mentioned process described in relation withFIG. 24 can be applied to the Home Price Index (“HPI”) as shown in FIGS.23A-23D to detect certain economic conditions. FIGS. 23A-23D show thescaled HPI (blue curve), the scaled lag-1 auto-correlation of AR(1)coefficients (green curve), and the scaled smoothed derivative of thelag-1 auto-correlation of AR1 coefficients (red curve). The results areshown for Charlotte in FIG. 23A, Las Vegas in FIG. 23B, Cleveland inFIG. 23C, and Phoenix in FIG. 23D. The last monthly calendar time beingused is August 2009. In this example, the scaled lag-1 auto-correlationof AR(1) coefficients and the scaled smoothed derivative of the lag-1auto-correlation of AR(1) coefficients can serve as the leadingindicators for HPI. A monitoring process can be configured to monitorthe leading indicators for the HPI.

In certain examples a computer-implemented method can be employed toperform one or more aspects of the methods discussed herein.Alternatively an apparatus comprising: a processor; and a memory forstoring computer readable instructions that, when executed by theprocessor, can cause the apparatus to perform one or more aspects of themethods discussed herein. In other embodiments, one or morenon-transitory computer-readable media may have instructions storedthereon that, when executed, cause at least one computing device toperform one or more aspects of the methods discussed herein.

Features of Methods and Systems for Determining Economic ConditionsAccording to Examples of the Disclosure:

In one example, a method can include one or more of the following steps:receiving an input of economic data including a plurality ofmacroeconomic and microeconomic indicators, applying a band pass filterto the economic data, extracting a cycle signal from the band passfilter, determining a synchronization event of the plurality ofmacroeconomic and microeconomic indicators, and outputting an economicprediction based on the synchronization event. The band pass filter canbe applied to the plurality of macroeconomic and microeconomicindicators simultaneously. After applying a band pass filter to theeconomic data the absolute value of the plurality of the filteredeconomic data series is calculated. The absolute value of the filtereddata can be reviewed to determine correlations and concentrationmeasures across the plurality of economic indicators to predict economictrends. A principal component analysis can be applied to the economicdata to determine a market driver.

In another example, a method can include one or more of the followingsteps:

determining a first principal component number from the economic data,determining eigenvector contents of the first principal componentnumber, normalizing the eigenvector contents of the first principalcomponent number, sorting the normalized eigenvector contents fromsmallest to largest, calculating the distances between the normalizedeigenvector contents, calculating the mean of the distances to obtain adispersion level, and outputting an economic instability predictionbased on the dispersion level of the distances. The predetermined timewindow can be set to approximately 70 months, though it can have anyvalue deemed reasonable. The first principal component number canrepresent systemic risk. The eigenvalue contents can be normalized suchthat the eigenvalue contents total summation is 100%. A lower dispersionlevel may indicate a greater risk of economic instability. The input ofeconomic data can include both microeconomic data and macroeconomicdata. Prior to applying principal component analysis the economic datacan be passed through a band pass filter to determine a cycle state.

In another example, a method can include one or more of the followingsteps: receiving economic data of a predetermined time window andapplying principal component analysis to the economic data, determininga first principal component number, determining a second principalcomponent number, calculating the ratio between the first principalcomponent number and the second principal component number, determiningthe largest ratio between the first principal component number and thesecond principal component number, and outputting a market driverprediction based on the largest ratio between the first principalcomponent number and the second principal component number. Thepredetermined time window can be approximately 70 months. The firstprincipal component can be systemic risk and the second principalcomponent can be non-systemic risk. The input of economic data mayinclude both microeconomic data and macroeconomic data.

In another example, a method can include one or more of the followingsteps: receiving economic data of a multiple time series, determiningcorrelation coefficient values using the data of multiple time series,creating a matrix using the correlation coefficient values, determiningthe nearest neighbor spacing distribution, determining lowest values ofcorrelation coefficients and set the lowest value of correlationcoefficients iteratively to zero, determining a threshold level by usinga chi-square test, and outputting an economic instability predictionbased on the threshold level. The nearest neighbor spacing distributioncan be determined by taking the differences between nearest neighboreigen-values derived from a correlation coefficient matrix. Thethreshold level can be determined by detecting a sharp phase transitionfrom a Gaussian Orthogonal Ensemble distribution to a Poissondistribution. A standard Pearson correlation coefficient values can bedetermined using the data of multiple time series.

In another example, a method can include one or more of the followingsteps: obtaining economic data of multiple time series of interest,determining correlation coefficient values using the data, forming amatrix of the correlation values and determining the eigen-values of thecorrelation matrix, forming a random matrix counterpart with knowneigen-value density distribution, comparing statistics of theeigen-values of the data correlation matrix to the values of the randommatrix counterpart of known eigen-value distribution, determining thedifference between the number of the eigen-values of the datacorrelation matrix and the random matrix counterpart, and outputting athe number of principal components to retain based on the differencebetween the number of the eigen-values of the data correlation matrixand the random matrix counterpart. The standard Pearson correlationcoefficients can be obtained using the data in the method.

In another example, a method can include one or more of the followingsteps: receiving data of a predetermined time series, determining autoregressive model coefficients for the data of the time series,determining a scaled lag-1 auto-correlation of the coefficients, andoutputting an indicator based on the scaled lag-1 auto-correlation ofthe coefficients. The indicator can be positive when the scaled lag-1auto-correlation of the coefficients is above a predetermined thresholdvalue. The indicator can be negative when the scaled lag-1auto-correlation of the coefficients is below a predetermined thresholdvalue. The data of the predetermined time series can be Home Price Indexdata. A scaled smoothed derivative from the scaled lag-1auto-correlation of the coefficients can be determined. A secondindicator can be output based on a scaled smoothed derivative from thescaled lag-1 auto-correlation of the coefficients.

As previously discussed, many complex dynamical systems have criticalthresholds called tipping points where the system abruptly shifts fromone state to another. In global finance, there is big concern about asudden systemic decline in values that can threaten the global economicsystem which, in turn, may lead to global economic downturns. It istraditionally difficult to predict such critical transitions since thestate of the system may show little change before the tipping point isreached. However, it seems that certain generic signs may take place ina wide class of systems as they approach a critical point. The dynamicsof systems near a critical point have generic properties, regardless ofdifferences in the details of each system. Critical thresholds for suchtransitions correspond to bifurcations. In the so-called catastrophicbifurcation, once a threshold is exceeded, a positive feedback mechanismpushes the system through a phase of directional change towards acontrasting state. The transitions from a stable equilibrium to a cyclicor chaotic state can also happen in some other types of bifurcations.Three possible early-warnings in the dynamics of a system approaching atipping point may be used in emerging instability detection: slowerrecovery from perturbations, increased autocorrelation and increasedvariance.

When an economic system approaches a critical threshold (e.g., a tippingpoint), the system may abruptly shift from one state to another (e.g.,where the value of underlying economic instrument becomes appreciablylower or higher than a previous value). Consequently, it may beimportant that an early warning signal be expeditiously conveyed to auser who is responsible for the underlying economic instrument. Those ofskill in the art will appreciate that there may be instances whereinstability may occur within minutes or seconds after receiving an earlywarning signal. Those of skill in the art will also appreciate thateconomic instruments may be a tradable resource of any kind. Economicinstruments may be cash, evidence of an ownership interest in an entity,a contractual right to receive or deliver cash or another economicinstrument, or a derivative. Cash instruments may be securities whichare readily transferable or may be instruments such as loans, where bothborrower and lender have to agree on a transfer. Derivative instrumentsmay derive their value from the value and characteristics of one or moreunderlying entities such as a resource, index, or interest rate.Derivative instruments may be exchange-traded derivatives andover-the-counter derivatives.

While the term “network” includes not only a “physical network” and a“content network,” as previously discussed, the term “network” may alsorefer to “platform networks.” For example, a Unified Network Platform(UNP) comprises an open platform architecture that enables theconvergence of network applications including data networking, security,network infrastructure, network applications, and network management.

FIG. 25 illustrates apparatus 2500 that comprises data processor 2501and early warning system 2502 in accordance with an aspect of theembodiments. For example, processor 2501 obtains information about areal-world system (which is typically external to apparatus 2500) suchas an underlying complex economic network. The information is gatheredfrom information source 2530, which may comprise a network of computersand/or databases that monitor the real-world system. As previouslydiscussed, processor 2501 processes the obtained information todetermine if the underlying real-world system exhibits characteristicsare indicative that the real-world system is approaching a criticalthreshold. In such a situation, early warning system 2502 generates anearly warning signal to one or more computing devices 2520-2522 toinform devices 2520-2522 about the potential instability of the system.

Those of skill in the art will appreciate that apparatus 2500 mayprocess and transform information from different types of real-worldsystems in the realms of biology, epidemiology, genetics,transportation, computer software, and economics.

Referring to FIGS. 23A-D, early warning system 2502 may determine fromthe scaled smoothed derivative that the home price index (HPI) for LasVegas (a geographical region as shown in FIG. 23B) is approaching acritical threshold while the HPI for Cleveland (a geographical region asshown in FIG. 23C) is not. In such a situation, early warning system2502 may generate an early warning signal having a positive warningsignal component corresponding to Las Vegas and a negative warningsignal component for Cleveland. For example, the early warning signalmay comprise a message having indicators for Las Vegas and Clevelandwith positive and negative values, respectively.

An early warning signal may comprise one or more messages that are sentto one or more computing devices. For example, as will be discussed, afirst message may notify computing devices 2520-2522 that a criticalthreshold is being approached while a second message has more detailedinformation about the possible instability. Based on the informationcontained in the early warning signal, action via a computing device maybe initiated to address the possible instability.

FIG. 26 illustrates computer system 2600 that supports apparatus 2500for processing information from information source 2630 and generatingan early warning signal pertaining to the stability/instability of areal-world system in accordance with an aspect of the embodiments.Computer system 2600 comprises servers 2601, 2602, and 2603. With someembodiments, data server 2601 obtains information about a real-worldsystem from information source 2630. As previously discussed, thereal-world system may comprise an underlying economic network thatexhibits structural attributes that may provide clues about economicinstability or economic stability. Based on processing the informationfrom information source 2630, data server 2601 may detect early-warningsigns of the real-world system approaching a critical threshold wherethe system becomes unstable. If so, data server 2601 generates an earlywarning signal. Notification of the event may be sent to one or moreusers (associated with computing devices 2620-2622) via communicationserver 2603.

Communication server 2603 may send an event notification to some or allof computing devices 2620-2622 through one or more differentcommunication media. For example, a notification may be sent via a shortmessage service (SMS) and/or via an intranet/internet. Referring to FIG.26, communication server 2603 interacts with primary computing device2620 via local communication network 2610 and with remote computingdevices 2621-2622 via external communication network 2611. With someembodiments, primary computing device 2620 may support a user when theuser is working in an office environment. However, the user may not beavailable in the office environment when the real-world system isexhibiting characteristics of possible instability and immediate actionmust be initiated by the user. In such a situation, the eventnotification may be sent to remote computing device 2621 and/or remotecomputing device 2622, which are registered to the user. As will bediscussed, if an acknowledgement to the notification is received bycommunication server 2602 from one of remote computing devices2621-2622, more detailed information about the possible instability maybe sent to the responding computing device. For example, ifcommunication server 2602 receives an acknowledgement from device 2622,an early warning application may be manually or automatically launchedat device 2622 in which a communication connection is establishedbetween web server 2603 and device 2622 via external communicationnetwork 2611. One or more messages may be sent to device 2622 via theearly warning application through the communication connection toprovide additional information about the detected event. For example,textual information identifying the corresponding geographical area(e.g., Charlotte, Las Vegas, Cleveland, and/or Phoenix as shown in FIGS.23A-D, respectively) and including quantitative information (e.g., ARauto-correlation and/or smoothed derivative of the AR) as well asgraphical information (e.g., plots showing relationships between varioustypes of data as shown in FIGS. 23A-D) may be sent from web server 2603to device 2622 via communication server 2602 and external communicationnetwork 2611. Detailed information may include information specific to aparticular geographical region associated with a registered device2620-2622 as well as information about other geographical regions thatmay be correlated with the particular geographical region. For example,a potential instability for Las Vegas may be detected so that detailedinformation may be sent a device registered for a user responsible forthe Las Vegas geographical region. However, the detailed information mayinclude HPI information about Charlotte, Cleveland, and/or Phoenixand/or other related information (e.g., home liquidity) so that theregistered user for Las Vegas region can correlate the potentialinstability with other geographical regions to ascertain whether thereis a systemic instability for the HPI.

With some embodiments, the detailed information in the early warningsignal may indicate whether the economy associated with the instabilityis moving in a positive or negative direction and further indicate ifthere are any particular elements driving the change. On a macro levelthis could be certain industries or asset classes. On a micro level itcould be a certain company's economic performance or group of companies'economic performances. Also, economic instability may be affected bypolitical instability, e.g., government grid-lock, elections, and thelike. As previously discussed, principal component analysis or PCA maybe used to identify the particular “drivers” of the trends. This canindicate which subcomponents are accounting for the change in state. Forexample, as previously discussed in reference to FIG. 9, largest ratiosof two predetermined components may illustrate a particular economicdriver. Two subcomponents can be compared by calculating the ratiobetween their first and second principal components. In one example, thefirst principal component can represent systemic risk, for example, arecession or other event external to the data being studied. The secondprincipal component can be non-systemic (idiosyncratic) risk, which canbe specific to a particular subcomponent. Subcomponent 5 may be themarket driver because the ratio of the first principal component to thesecond principal component is the largest.

With some embodiments, an early warning application may be automaticallylaunched at computing device 2620, 2621, or 2622 whenever anacknowledgment is generated at the device in response to receiving anearly warning notification. Once launched, the early warning applicationmay extract additional information about a potential instability fromone or more messages contained in an early warning signal, appropriatelyfilter the additional information, format the processed information, anddisplay the information at the computing device. Moreover, the earlywarning application may be device-independent so that the applicationcan execute on a variety of different device types.

Servers 2601, 2602, and 2603 may be combined on the same physicalmachines, and retain separate virtual or logical addresses, or mayreside on separate physical machines. FIG. 26 illustrates just oneexample of a network architecture that may be used, and those of skillin the art will appreciate that the specific network architecture anddata processing devices used may vary, and are secondary to thefunctionality that they provide, as further described herein. Forexample, services provided by data server 2601, web server 2602, andcommunication server 2603 may be supported by a single physical server.

Each computing component (servers 2601-2603 and computing devices2620-2622) may be any type of known computer, server, or data processingdevice as discussed herein, e.g., as shown in FIG. 1.

FIG. 27 illustrates flow diagram 2700 for generating an early warningsignal from processed information and responding to the early warningsignal in accordance with an aspect of the embodiments. With someembodiments, flow diagram 2700 may be performed by one or moreprocessors executing computer readable instructions from one or morenon-transitory computer-readable media spanning data server 2601,communication server 2602, and/or web server 2603.

Referring to FIGS. 26 and 27, system 2600 obtains information about areal-world system from information source 2630. At block 2701, system2600 determines whether there is a potential instability within thereal-world system. For example, system 2600 may detect the dynamics of asystem approaching a tipping point as indicated by slower recovery fromperturbations, increased autocorrelation, and/or increased variance. Ifnot, system 2600 continues to monitor/analyze information from source2630.

If a potential instability is detected at block 2701, system 2600determines the assigned registered devices for the potential instabilityat block 2703, where registered devices may be associated with one ormore users. For example, the real-world system may exhibit instabilityat some or all of the associated geographical regions (entities), whereregistered devices may be assigned to a particular geographical region.Based on the determined geographical region, one or more registeredcomputing devices 2620-2622 (e.g., desktops, laptops, notebooks, mobiletelephones or smartphones with applications and other functionality,handheld device with Wi-Fi or other wireless enabled tablets, tabletcomputers, PDAs, and the like), may be registered for the geographicalregion on behalf of the user or users responsible for the identifiedgeographical region.

Devices 2620-2622 may be registered according to different criteria,e.g., associated with specific geographical regions as discussed above.Moreover, embodiments may consider an identified economic systemspanning geographies, industries, and/or other economic neighborhoodswhen determining criteria for the registration of devices 2620-2622.

System 2600 determines whether any of registered computing devices2620-2622 are active at block 2704. If so, system 2600 generates anearly warning notification to the primary device at block 2705. Forexample, a user may be active at primary computing device 2620 (whichmay be registered as the primary device since it operates in a moresecure office environment). With some embodiments, the notification maybe a short text message identifying the possible instability and theassociated geographical region.

At block 2706, an early warning application may be automaticallylaunched at the primary device after the notification is received at thecomputing device. The early warning application may obtain detailedinformation about the detected event. For example, a communicationconnection may be established between web server 2603 and device 2620via communication server 2602 to provide quantitative information (e.g.,AR auto-correlation and/or smoothed derivative of the AR) as well asgraphical information. The additional information may be formatted inone or more webpages by web server 2603 and displayed at computingdevice 2620 via the early warning application and/or web browser.

However, if system 2600 determines that none of the registered devicesis active at block 2704, system 2600 generates an early warningnotification to all registered devices at block 2707. For example, ifdevice 2620 is not active, system 2600 may send an early warningnotification to remote computing devices 2621-2622, which may beportable devices having lesser capabilities than device 2620.

If system 2600 receives an acknowledgement from a registered device, asdetected at block 2708, the early warning application may beautomatically launched at the responding device at block 2709. Inaddition, based on the display capabilities of the device, some or allof the currently displayed framed graphical elements (sometimes referredas windows) may be modified, resized, and/or minimized to moreeffectively display the detailed content of the early warning signal.For example, a currently displayed window may be obscuring a windowdepicting a graphical representation about the HPI (FIGS. 23A-D). If so,the currently displayed window may be modified or moved to reduce theobscurity.

If system 2600 does not receive an acknowledgement to the early warningnotification from any of the registered devices, system 2600 mayinitiate an alternative early warning notification. A notification maybe sent to alternative devices (e.g., devices registered for a differentgeographical region) or via e-mail sent to registered users at block2710.

FIG. 28 illustrates flow diagram 2800 for generating an early warningsignal from processed information and responding to the early warningsignal in accordance with an aspect of the embodiments. Referring toFIG. 27, blocks 2801-2808 and 2810 are similar to blocks 2701-2708 and2710, respectively. However, if system 2600 receives an acknowledgmentfrom a registered remote device at block 2808, the early warningapplication is launched at the primary computing device (e.g., theuser's computer located on the office environment) to retrieve detailedinformation from web server 2603. A communication connection may then beautomatically established between the remote device and the primarydevice via communication server 2602. Consequently, a user may obtaindetailed information about the potential instability from the primarydevice via the communication connection at the remote device.

Referring to FIGS. 27 and 28, while early warning notifications maytarget specific devices, embodiments may target access-id applications,e.g., social media user accounts, domain-specific registered users, andthe like. The additional layer of notification may not have veryspecific details of the alert other than notify the recipient of themessage to check the devices where the recipient chooses to receive thesecured alert.

FIG. 29 illustrates flow diagram 2900 for sending an early warningsignal to a plurality of computing devices in accordance with an aspectof the embodiments. As with flow diagrams 2700 and 2800, system 2600generates an early warning notification to all registered computingdevices 2620-2622 at block 2901. When the first device acknowledges thenotification at block 2902, system 2600 provides detailed early warninginformation to the first responding device. At block 2903, an invitemessage may be generated and sent to another device so that the user ofthe first device can collaborate with another user of the other device,where the other device may or may not be registered. When communicationis established between the first device and the other device, thedetailed information may be shared with the other device. Furthermore, asubsequent invite massage may be generated by the first device, otherdevice, or system 2600 to include a further device during thecollaboration.

FIG. 30 illustrates flow diagram 3000 for displaying information derivedfrom an early warning signal at a computing device in accordance with anaspect of the embodiments. At block 3001, system 2600 receives anacknowledgment to an early warning notification from one of theregistered devices. System 2600 identifies the device type of theresponding device at block 3002 and determines the displaycharacteristics based on the identified device type at blocks 3003-3005,which corresponds to a desktop, smartphone, and tablet, respectively.However, additional or different types of devices may be supported.Based on the obtained display characteristics, web server 2603 mayconsider the viewable image size, resolution, and aspect ratio todetermine the amount of detailed information provided to the device. Webserver 2603 may partition the detailed information into a plurality ofinformation components and provide selected information components ofthe detailed information to the device based on the displaycharacteristics at block 3006. For example, different informationcomponents may comprise different graphical elements having varyingdegrees of information and/or resolution.

Information components of the detailed may span a variety of formsincluding textual information and visualizations that may convey thesystematic uncertainty associated with the potential instability.Visualizations may be presented in different forms, including graphs,network diagrams, correlograms, and the like. The information componentsmay be pushed from computer system 2600 or may be pulled by computingdevices 2620-2622. When pushed, computer system 2600 may obtain thedisplay characteristics of a responding computing device 2620-2622, andsend information components that the responding computing device caneffectively display. When pulled, an early warning application maysupport a dashboard from which a user can select available informationcomponents that are to be displayed on the computing device.

When early warning system 2502 has notified registered devices about thepotential instability, early warning system 2502 may further insure thatthe computing devices of appropriate individuals are connected via acommunication system (for example, by generating invite messages asdiscussed above) so that a decision about the instability can be reachedby a participant or collectively by a group of participants. With someembodiments, a data structure (e.g., role-based matrix) may map a levelof responsibility for each of the participants and with identities ofassociated computing devices for those who are needed in thedecision-making process. The data structure may be based on anorganization chart, where typically the higher one's level inorganization chart, the greater the authority one has for making adecision. The data structure may relate the individual, the level ofresponsibility, and the device identity of the individual. The datastructure may reside within system 2600 or may reside in a separatecomputing entity (e.g., a separate server).

With some embodiments, early warning system 2502 may interact with anemergency notification system that may have been established to reachpersonnel in the event of natural and non-natural events. The emergencynotification system may follow a contact tree to contact employees viatelephone calls, e-mails, and so forth. The early warning signal fromearly warning system 2502 may induce the early notification system tofunction as an added contact mechanism. Moreover, the contact mechanismmay include delivering early warning signals through social media andsoftware-defined networks.

In some situations, only some of the participants of the communicationconnection may have the authority to make a decision about theinstability. With some embodiments a separate communication connectionmay be established among the decision-makers (e.g., a subset of theparticipants) so that messaging/conversation among the decision-makersis not available to other participants. For example, some of theparticipants in a first communication connect may not have the authorityto make a decision about an instability. If so, a second communicationconnection may be established among the decision-makers, where theconversation/messaging of the second communication connection is notavailable to the first communication connection.

Although the subject matter has been described in language specific tostructural features and/or methodological acts, it is to be understoodthat the subject matter defined in the appended claims is notnecessarily limited to the specific features or acts described above.Rather, the specific features and acts described above are disclosed asexample forms of implementing the claims.

What is claimed is:
 1. An apparatus comprising: an early warning system;a processor; and memory storing computer readable instructions that,when executed by the processor, cause the apparatus to: receive economicdata of multiple time series for an economic system; determinecorrelation coefficient values using the data of the multiple timeseries; create a matrix using the correlation coefficient values;determine the nearest neighbor spacing distribution; determine lowestvalues of correlation coefficients and set the lowest values ofcorrelation coefficients iteratively to zero; determine a thresholdlevel of a control parameter, wherein the threshold level is detectedonce a phase transition occurs and wherein a modeled distribution forthe economic data changes from a first distribution to a seconddistribution; output an economic instability prediction based on thethreshold level; generate, based on the economic instability prediction,an early warning signal about the economic data, wherein the earlywarning signal is indicative of whether the threshold level isapproaching and wherein the early warning signal comprises an earlywarning notification and detailed information about a potentialinstability of the economic data; determine a set of registeredcomputing devices associated with the economic system, wherein the setincludes a first device and a second device; send the early warningnotification about the economic system to the first and second devices;initiate launching of an early warning application at the first device;and send the detailed information about the potential instability to thefirst device via the early warning application.
 2. The apparatus ofclaim 1 wherein the nearest neighbor spacing distribution is determinedby taking the differences between nearest neighbor eigen-values derivedfrom a correlation coefficient matrix.
 3. The apparatus of claim 2,wherein the memory storing the computer readable instructions that, whenexecuted by the processor, further cause the apparatus to: diagonalizethe correlation coefficient matrix until the phase transition occurs. 4.The apparatus of claim 1 wherein the threshold level is determined byusing a chi-square test.
 5. The apparatus of claim 4 wherein thethreshold level is determined additionally by detecting a sharp phasetransition from a Gaussian Orthogonal Ensemble distribution to a Poissondistribution.
 6. The apparatus of claim 1 wherein standard Pearsoncorrelation coefficient values are determined using the data of multipletime series.
 7. The apparatus of claim 1, wherein the memory storescomputer readable instructions that, when executed by the electronicprocessing device, further cause the apparatus to: send an invitemessage to another device, wherein the invite message invites the otherdevice to communicate with the first device.
 8. The apparatus of claim1, wherein the memory stores computer readable instructions that, whenexecuted by the electronic processing device, further cause theapparatus to: when the acknowledgment to the warning notification isreceived from the first device, identify a device type of the firstdevice; determine a first display characteristic of the first devicefrom the device type; and send a subset of information components to thefirst device based on the first display characteristic, wherein thedetailed information comprises a plurality of information components. 9.The apparatus of claim 1 further comprises: a data structure, whereinthe data structure maps decision-making roles of individuals withassociated device identities, wherein the device identities identify aplurality of devices; and wherein the memory device storing computerreadable instructions that, when executed by the electronic processingdevice, further causes the apparatus to establish a first communicationconnection with the plurality of devices based on the data structure.10. The apparatus of claim 9 wherein the memory stores computer readableinstructions that, when executed by the electronic processing device,further cause the apparatus to: establish a second communicationconnection with a subset of the plurality of devices based on the datastructure; and prevent content of the second communication connectionfrom being conveyed to the first communication connection.
 11. Theapparatus of claim 1 wherein the memory stores computer readableinstructions that, when executed by the electronic processing device,further cause the apparatus to: send an alert indication to a designatedrecipient through an access-id application, wherein the recipient isinformed to check at least one of the set of registered computingdevices.
 12. A method comprising: receiving, by a processing system,economic data of multiple time series; determining, by the processingsystem, correlation coefficient values using the data of the multipletime series; creating, by the processing system, a matrix using thecorrelation coefficient values; determining, by the processing system,the nearest neighbor spacing distribution; determining, by theprocessing system, lowest values of correlation coefficients and settingthe lowest values of correlation coefficients iteratively to zero;determining, by the processing system, a threshold level of a controlparameter, wherein the threshold level is detected once a phasetransition occurs and wherein a modeled distribution for the economicdata changes from a first distribution to a second distribution;outputting, by the processing system, an economic instability predictionbased on the threshold level; generating, by the processing system, anearly warning signal based on the economic instability prediction for aneconomic system, wherein the early warning signal is indicative ofwhether the threshold level is approaching and wherein the early warningsignal comprises an early warning notification and detailed informationabout a potential instability of the economic system; sending the earlywarning notification about the economic system to a primary device; whenan acknowledgement is received from the primary device, sending thedetailed information about the instability to the primary device; andwhen the acknowledgement is not received from the primary device,determining a set of registered computing devices associated with theeconomic system, wherein the set includes a first device and a seconddevice; sending the early warning notification about the economic systemto the first and second devices; initiating a launching of an earlywarning application at the primary device; remotely connecting theprimary device to the first device; and sending the detailed informationto the first device via the primary device.
 13. The method of claim 12wherein the nearest neighbor spacing distribution is determined bytaking the differences between nearest neighbor eigen-values derivedfrom a correlation coefficient matrix.
 14. The method of claim 12wherein the threshold level is determined by using a chi-square test.15. The method of claim 14 wherein the threshold level is determined bydetecting a sharp phase transition from a Gaussian Orthogonal Ensembledistribution to a Poisson distribution.
 16. The method of claim 12wherein standard Pearson correlation coefficient values are determinedusing the data of multiple time series.
 17. One or more non-transitorycomputer-readable media having instructions stored thereon that, whenexecuted, cause at least one computing device to: receive economic dataof multiple time series for an economic system; determine correlationcoefficient values using the data of the multiple time series; create amatrix using the correlation coefficient values; determine the nearestneighbor spacing distribution; determine lowest values of correlationcoefficients and set the lowest values of correlation coefficientsiteratively to zero; determine a threshold level of a control parameterby using a chi-square test, wherein the threshold occurs at an observedphase transition and wherein a modeled distribution for the economicdata changes from a first distribution to a second distribution; outputan economic instability prediction based on the threshold level;generate an early warning signal from the economic instabilityprediction for the economic system, wherein the early warning signal isindicative of whether the threshold level is approaching and wherein theearly warning signal comprises and early warning notification anddetailed information about a potential instability of the economicsystem; determine a set of registered computing devices associated withthe economic system, wherein the set comprises a first device and asecond device; send the early warning notification about the economicsystem to the first and second devices; initiate of an early warningapplication at the first device; and send the detailed information aboutthe potential instability to the first device via the early warningapplication.
 18. The one or more non-transitory computer-readable mediaof claim 17 wherein the nearest neighbor spacing distribution isdetermined by taking the differences between nearest neighboreigen-values derived from a correlation coefficient matrix.
 19. The oneor more non-transitory computer-readable media of claim 17 wherein thethreshold level is determined by detecting a sharp phase transition froma Gaussian Orthogonal Ensemble distribution to a Poisson distribution.20. The one or more non-transitory computer-readable media of claim 17wherein standard Pearson correlation coefficient values are determinedusing the data of multiple time series.
 21. An apparatus comprising: aprocessor; and memory storing computer readable instructions that, whenexecuted by the processor, cause the apparatus to: obtain economic dataof multiple time series of interest about an economic system; determinecorrelation coefficient values using the data; form a matrix of thecorrelation values and determine eigen-values of the correlation matrix;form a random matrix counterpart with a known eigen-value densitydistribution; compare statistics of the eigen-values of the datacorrelation matrix to values of the random matrix counterpart of theknown eigen-value distribution; determine a difference between aquantity of the eigen-values of the data correlation matrix and therandom matrix counterpart; output a retained number of principalcomponents to be retained based on the difference between the quantityof the eigen-values of the data correlation matrix and the random matrixcounterpart; and determine, from the retained number, a driver of trendsdetermined from the economic data; include an indicator about the driverin an early warning signal, wherein the early warning signal comprisesan early warning notification and detailed information about a potentialinstability of the economic system; determine a set of registeredcomputing devices associated with the economic system, wherein the setincludes a first device and a second device; initiate a launching of anearly warning application at the first device; and send the detailedinformation about the potential instability to the first device via theearly warning application.
 22. The apparatus of claim 21 whereinstandard Pearson correlation coefficients are obtained using the data.23. The apparatus of claim 21, wherein the memory storing the computerreadable instructions that, when executed by the processor, furthercause the apparatus to: determine maximum values of observedeigen-values; determine an observed quantity of observed eigen-valuesgreater than the maximum values; and approximate the retained number ofprinciple components by the observed quantity.
 24. The apparatus ofclaim 21, wherein the memory storing the computer readable instructionsthat, when executed by the processor, further cause the apparatus to:analyze micro level data after the driver has been determined.